Anisotropy in spatial order-disorder transformations and the 
vortex lattice symmetry transition in YNi2B2C and 

LuNi2B2C 

D. Jaiswal-Nagar^'^'^*, T. Isshiki^, N. Kimura^, H. Aoki^, H. Takeya"^, S. Ramakrishnan^ 
and A.K. Grover^ 

^Centre for Low Temperature Science, Tohoku University, Sendai, Miyagi, 980-8578, Japan 
^DCMP&MS, Tata Institute of Fundamental Research, Mumbai-400005, India 
^Physikalishes Institut, Goethe- Universitat, 60438 Frankfurt am Main, Germany 
■^National Institute for Materials Science, Sengen 1-2-1, Tsukuba, Ibaraki 305-0047, Japan 
* e-mail : Jaiswal-Nagar @physik. uni-fr ankfurt . de 

Abstract 

Explorations of the order-disorder transformation in vortex matter in single crystals of 
tetragonal structured (c/a ~ 3) borocarbide superconductors, YN 12620 and LuNi2B2C , reveal 
that vortex arrays experience different effective pinning in different crystallographic directions. 
We surmise that correlation exists between the large anisotropy in effective pinning/disorder 
and the differences in the (local) symmetry transition from rhombohedral to (quasi) square 
vortex lattice(VL). For field along high symmetry directions, like, c-axis and ab-basal plane, 
the VL symmetry is close to square and the ordered state spans a large field interval. When 
the field is turned away from the c-axis towards ab-plane, at intermediate angles, the region of 
ordered state shrinks, in response to enhancement in effective pinning. At such intermediate 
angles the symmetry of the VL would be far from ideal triangular or square. 

Introduction 

Disorder can play an instructive role in unraveling important aspects of phase transitions in 
a condensed matter system. In the context of vortex arrays, the role of two kinds of dis- 
order, namely, thermal fluctuations and quenched random disorder (QRD), in enunciating 
different phases of vortex matter is well documented [H [21 [3l IH O [6]. Quaternary boro- 
carbide superconductors have emerged as convenient test beds to explore a variety of newer 
issues in vortex physics, e.g., field induced changes in the symmetry of the vortex lattice 
(VL)[7l El m [ini ttH llSl [HI [Mj, symmetry changes due to interplay between non-local ef- 
fects in electrodynamics and symmetry of the superconducting order parameter |16i 1171 [TBI I19j. 
etc. From angle dependent dc magnetisation hysteresis and ac shielding response measurements 
in crystals of two borocarbides, viz., Y {Lu)N 12620 (Y(L)NBC), we show that the effective 
disorder experienced by the VL is different in different directions, which refiects in the details 
of their vortex phase diagrams. In particular, we conjecture that the closer the VL conforms to 
the perfect square/triangle symmetry, larger is the correlation volume for the ordered VL and 
it spans a larger portion of vortex phase space. 

Our study is based on a premise that a transformation/transition amongst vortex phases can, 
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prima facie, be captured via an anomalous change in the pinning attribute of the vortices, 
as evidenced in the field-temperature (H,T) dependence of the critical current density Jc of 
a given superconductor. We obtain information on Jc(H,T) via an exploration of the peak 
effect (PE) phenomenon (typically occurring near /prior to the normal-superconductor phase 
boundary [SIHIE]) and a second magnetization peak (SMP) anomaly [316] lying deeper in the 
mixed state of weakly pinned single crystals of YNBC and LNBC. Both these compounds are 
prime examples of field induced transitions in the local symmetry of the VL. For H || c, VL 
in them [HI [9l \TU[ \TT\ \T2\ [T3l I14j . initially undergoes a first order transition from a rhombus 
symmetry of one kind {(3 < 60°) to another (/3 > 60°) via a 45° reorientation, which eventu- 
ally becomes a square symmetry via a second order transition with progressive increase in field 
[H m [lOl [TTl I12j . For H _L c, however, the first order transition is via a 90° reorientation and 
the rhombus (/3 > 60°) never reaches a square symmetry (saturates with /? ~ 82°) [13\ I14j . 
From angle dependent PE measurements, we find that the details of the symmetry transitions 
in the VL in different orientations result in deciding the quality of spatial order in VL. The 
orientation of the magnetic field in directions along which the VL probably does not assume 
a symmetric shape (triangle/square with /3 — > 60°/90°, respectively), the vortex state is only 
partially ordered, and it eventually amorphizes in multiple steps via the occurrence of SMP 
and/or multiple magnetization peaks (MMP) apart from the PE. To fortify our assertion that 
the anomalies prior to the PE are (purely) disorder induced, we enhance the driving effect of 
the ac field on the partially ordered VL to improve its quality, which in turn suppresses the 
occurrence of anomalies precursor to PE. 

Experimental details 

The YNBC crystal (Tc ~ 15 K), grown by using an infrared image furnace [15], was cut as 
parallelopiped (see Fig. 1(b)) with dimension 1 = 3 mm, b = 0.7 mm, t = 0.67 mm. Its largest 
dimension is parallel to the a-axis of YNBC. The LNBC crystal on the other hand, was grown by 
the flux method [3; it is in the form of a thin platelet, with the c-axis perpendicular to its plane. 
The ac susceptibility measurements were performed on a top loading dilution refrigerator at the 
Centre for Low Temperature Science, Tohoku University, with amplitudes of the ac field ranging 
from 0.5 mT to 10 mT (peak to peak). Before measuring the ac susceptibility response, the 
sample was aligned along a particular crystallographic direction with the help of Laue diffraction 
method, within an accuracy of it 2°. For angular dependent studies, the crystal was rotated (i) 
from [100] to [110] in the (001) plane, (ii) from [110] to [001] and vice- versa in the (110) plane and 
(iii) from [001] to [100] and vice- versa in the (010) plane. It was noted that observed response at 
a given angle, was independent of the phase of the rotation. DC magnetization measurements 
were made using a commercial 12 Tesla Vibrating Sample Magnetometer (Oxford Instruments, 
U.K.) at the Tata Institute of Fundamental Research. The upper critical fields, Hc2, for LNBC 
crystal at T ~ 0.1 K, for H || [001], [110] and [100] are estimated to be 7.4 T, 7.9 T and 8.7 T, 
respectively. Each of these values corresponds to an ac susceptibility response merging into 
background. 

Results and discussion 

The inset panel of Fig. 1 (a) shows a M-H hysteresis loop in YNBC at 2.06 K with H [| c. The 
main panel of Fig. 1 (a) depicts the data on an expanded scale (blue curve for H [ | c and black 
curve for H _L c ) to focus on the PE anomaly seen as a hysteresis "bubble". It can be noted 
that while the magnetisation hysteresis widths (oc AM (H)) are notionally comparable in the 
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disordered PE region for both the orientations, they are very different in the field range prior 
to the onset of the PE (H < Hon)', is neghgibly small for H || c as compared to that for 
H _L c. Recall that AM (H) cx Jc (H) [201 121j , which in turn oc ^ [22], where Vc denotes the 
correlation volume over which the VL is well correlated. The above observation then, implies 
that a better ordered VL (with a large Vc) is formed when H 1 1 c, as compared to that for H ± c. 
This finding is significant as the superconducting parameters of tetragonal structured (c/a ~ 3) 
Y(L)NBC show little anisotropy (1 ~ 1.3) p3l[M[25]. 
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Figure 1: (Color online) (a) Main panel depicts a portion of the M-H hysteresis loop in a paral- 
lelopiped crystal of YNi2B2C for H ± c (black) and for H || c (blue), while the inset shows the full 
M-H loop for H || c. The onset field ifo„, the peak field Hp of the peak effect, as well as the upper 
critical field (at which the superconductor becomes normal) ifc2, have been marked for H _L c. (b) 
Normalized hysteresis width (oc Jc) vs. reduced field h = fioH / fiQHc2, corresponding to the M-H 
plots in panel (a). The dotted line near zero jc is a guide to the eye, it represents the experimental 
accuracy in the measured hysteresis width. 



During exploration of the PE phenomenon in spherical as well as an irregular shaped single 
crystals of a cubic superconductor VsSi, Kiipfer et al. [20] found that, prima facie, the disordered 
state of the VL had no angular dependence for H > Hp, while the state of order before the onset 
of PE crucially depends on the symmetry of the field direction with respect to the underlying 
CL. Taking a cue from this observation, we show in Fig. 1 (b), the plots of AM normalized to 
its value at peak field Hp, versus reduced field, h (= fioK/ fiQHc2) for H || c and H _L c. It is 
satisfying to note that both the curves merge for H > Hp (similar to those in the VsSi case [26j), 
thereby elucidating the insensitivity of the disordered vortex matter to the field-direction, for 
H > Hp, in YNBC as well. 

Considering that the superconducting parameters (coherence length id^ab, penetration depth 
^c/^ab, etc.) show little anisotropy (~ 1-1.2) for rotation from the c-axis to the basal ab-plane 

{AMh /AA/Hon) 

directions |23l[2ll[25], the observation that t-t-t-t — t^tt ~ 17, suggests that the behaviour 

' ' {A'VlHp/l^MHon) H^c 

of the spatial order of the VL is not governed by the intrinsic superconducting anisotropy of 
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YNBC. It is fruitful at this juncture to recall reports [HI [131 [H] that a square lattice gets formed 
for H II c, whereas only a quasi-square lattice (apex angle ~ 82°) is observed for H _L c. This 
observation along with the inference of a larger correlation volume for ordered VL for H || c 
compared to that for H _L c, then implies that a better spatially ordered VL is obtained, when 
its underlying symmetry is a perfect square rather than a quasi-square. 

In order to obtain information on the spatial ordering in the VL in directions other than H || 
c and H _L c, we took a recourse to ac susceptibility measurements in LNBC, on which dHvA 
measurements were already being performed at 0.1 K [27J. In the main panel of Fig. 2 we show 
angle dependent plots of normalised in-phase susceptibility Xat (= x' i^) / x' (Hp)) at T = 0.1 K 
with an hac of 1 mT (peak to peak) at 67 Hz. The normalisation for this as well as subsequent 
measurements has been chosen in accordance with the discussion above. The PE feature is 
clearly observable for all the angles in Fig. 2. Note first that prior to the PE region (H < Hon), 
Xjv (c< Jc) [29] is the smallest for H || c (i.e., [001]) as compared to the other orientations, 
suggesting that the VL of LNBC also has better spatial order for H 1 1 c as compared to any of 
the directions within the ab-basal plane. Also, note that within the basal plane, for H < Hon, 
Xjv is least for H || [110] and highest for H || [100], thereby indicating that the VL for H || [110] 
is better spatially ordered than that for H || [100], reminiscent of the behaviour reported in VsSi 
[26]. 

The inset panel in Fig. 2 shows plots of x'n reduced field h for rotation of the field from 
[100] to [110]. In terms of h, the onset field of the PE (hon), increases from ~ 0.65 for H || [100] 
to a value of ~ 0.72 for H || [110]. A premise, that hon is governed by the residual disorder in 
the ordered state of the VL, (i.e., higher the residual disorder, lower is the hon and vice-versa) 
[2H] . then, suggests that the vortex state for H || [110] is not only better ordered as compared 
to any other orientation in the ab-basal plane, but the ordered state also spans a larger field 
interval prior to the PE. This finding is significant as the intrinsic superconducting parameters 
of the tetragonal structured (c/a ~ 3) Y(L)NBC compounds show little anisotropy (1 ~ 1.2) 

[231 [HIES]. 

We now recall that Park et al. |16j had surmised that hexagonal to (quasi) square transition 
for VL in LNBC occurs at a lower threshold field for H || [110] as compared to that for H || [100]. 
This could imply that well ordered Bragg glass (BG) like phase [2] with no dislocations and 
having (quasi) square symmetry could permeate the entire sample at a lower field for H || [110] as 
compared to that for H || [100]. The occurrence of symmetry re-orientation transition at a lower 
field value facilitates the extension of spatial correlations to a larger volume. The symmetry 
transition in a VL is understood to arise due to a correlation of the VL with the underlying 
CL [17]. A non-local relationship between the super-current and the vector potential results 
in a distortion of the super-current, which leads to an anisotropic intervortex interaction. The 
threshold field at which non-local relationship governs the symmetry transitions in VL for H || 
c has been shown to enhance as the electron mean free path due to impurity effect progressively 
decreases in YiXi{Nii-xCox)2B2C series of alloys [18]. If such a notion is extended to hold for 
the state of order of the VL for H ± c, it would in turn confirm that the effective disorder being 
experienced for H || [110] is less than that for H || [100] in LNBC. 

To elaborate more on anisotropy in effective disorder, we show in Fig. 3 plots of Xat ^s. 
applied field, for its movement out of the ab-basal plane in the LNBC crystal. Panels (a)/(b) 
show seven plots each for a change in orientation of the applied field from [110]/ [100] (red curve) 
towards [001] direction (black curve), in steps of 15 degrees. The color coding and numerical 
labeling (1 to 7) identify the corresponding angles (measured w.r.t. [110]/[100]). Presence of 
a PE anomaly can be noted for all the orientations. As the field stands changed by 60° from 
[110] towards [001] (cf. violet curve - 5 in the panel (a)), one can witness multiple anomalies, 
termed as multiple magnetisation peaks (MMP), preceding the arrival of PE. However, such 
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Figure 2: (Color online) The solid curves in the main panel show the angular dependence of nor- 
malised in-phase ac susceptibility in LuNi2B2C at T = 0.1 K for the in plane rotation of the field 
from [100] to [110] in increments of 15°, while the black dotted curve represents the x'n H || [001]. 
The onset field of PE Hon-, the peak field Hp, the irreversibility field Hirr-, and the upper critical field 
Hc2 have been marked for H || [100] (blue curve). The inset shows x'n ^s. h (= fioB./ fioHc2) for the 
rotations in the basal plane from [100] to [110]. 
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Figure 3: (Color online) Plots of normalised in-phase ac susceptibility {x'n) applied field in 
LuNi2B2C at T = 0.1 K, for change in orientation of the applied field from (a) [110] to [001] and 
(b) [100] to [001], in increments of 15°. 



anomalies are not evident for the corresponding rotation from [100] in the panel (b) (cf. violet 
curve - 5 in panel (b)). On changing the field-angle to 45° in Fig. 3 (a)(cf. green curve - 4), the 
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PE anomaly remains considerably broadened and it can be imagined to be a juxtaposition of 
a SMP like anomaly preceding the PE. However, only a broad PE anomaly is observed for the 
corresponding rotation in Fig. 3 (b) (cf. green curves - 4). As the orientation of the applied field 
stands changed by 30° from [110]/ [100] (cf. pink curves - 3), only the PE anomaly is evident 
in the panel (a), but a SMP like anomaly now surfaces up in the panel (b). For the applied 
field changed by 15° from [110]/[100] (cf. blue curves - 2), the SMP anomaly can be imagined 
to move closer and juxtapose with PE in the panel (b), whereas only a PE is witnessed in the 
panel (a). For H || [110]/[100] (cf. red curves - 1), only a PE anomaly is observed, which as 
anticipated, is sharper in panel (a) vis-a-vis panel (b). 

SANS measurements on YNBC had revealed [9], that when the field was tilted 30° away from 
the c-axis along the ac plane, a distorted hexagonal VL {(3 ~ 75°) remained stabilised upto the 
highest measured field of 0.24 Tesla. In Fig. 3 (a), MMP anomalies are observed for such a 30° 
change from [001] direction (cf. violet curve - 5 in the panel (a)). Correspondence between SANS 
data in YNBC and our observation in Fig. 3 in LNBC could imply that when the applied field 
lies in a crystallographic direction that the symmetry of the VL is not close to (quasi) hexagonal 
with /? ~ 60° or (quasi) square with /? ~ 90°, the spatial correlations in the VL could stand 
compromised as if there exists larger effective disorder. For field-angles such as 45° w.r.t. [110] 
(curve - 4 in Fig. 3 (a)) or 30° w.r.t. [110] (curve - 3 in Fig. 3 (b)), there are no SANS data 
in the literature, however, we anticipate that the VL symmetry, that results in SMP anomaly 
preceding the PE anomaly, is not high in such directions. 
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Figure 4: Plots of (negative) in-phase ac susceptibility (-%') vs. applied field in LuNi2B2C at T = 
0.1 K, for H II 45° tilted from [110] towards [001], using an ac amplitude of 0.5 mT, 1 mT and 10 mT. 
The onset field (Hon), peak field {Hp), irreversibility field {Hirr) and upper critical field (-ffc2) have 
been marked for hac = 10 mT. 



To demonstrate that SMP anomaly is, (purely) disorder induced, the disordered VL was 
shaken with a higher ac field to drive out the precursor SMP anomaly. Fig. 4 depicts the ac 
susceptibility data (-%' vs. field for different amplitudes, hac, of the ac field) at T = 0.1 K for 
H corresponding to curve - 4 in Fig. 3 (a). Note that for hac = 0.5 mT (dashed line), and 
hac = 1 mT (solid line), both SMP anomaly and PE anomaly are present. However, when 
the applied ac field is raised to 10 mT, only the PE anomaly can be witnessed. The ac field 
not only invokes the shielding currents in a superconductor, but also acts as a driving force, 
by attempting to make the metastable vortices reach their equilibrium configuration, while 
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shaking them. The suppression of the SMP anomaly by shaking the VL estabhshes that the 
SMP anomaly is nucleated by an enhanced effective disorder. 
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Figure 5: The two panels show the variation of Hc2, Hp and Hon, when the applied field is changed 
from (a) [110] to [001] and (b) [100] to [001]. Dash and Dots are the guides to the eye, while the 
solid line is the fitting of Hc2{0) to the Ginzburg-Landau theory. Various phases, e.g., Bragg glass, 
polycrystalline vortex glass and amorphous vortex state have been named, as per nomenclature in 
the literature [21 [30]. 



Collating field-angles at which some typical (like, Hc2) as well as anomalous features (like, 
Hon and Hp) could be noted in the ac susceptibility, we show in Fig. 5 the vortex phase dia- 
grams for LNBC. The panels (a)/(b) in Fig. 5 correspond to the field tilted from [110]/[100] 
towards [001]. Dashes and dots are only guides to the eye, while the solid line is the fit to 
Hc2i0) using the Ginzburg-Landau relationship Hc2 = Hc2i\\c,T){sin'^ (6) + e'^cos'^ie))-'^/'^ , 
where e = -^^jp^y, and 9 is the angle between the applied field H and the ab-basal plane. 
From the fitting, we can obtain, e ~ 0.86 and 0.91, for field tilted from [110] and [100], respec- 
tively to [001], consistent with the earlier estimates [21] ■ The average out of plane anisotropy, 
0.5(i?|.2'^'^^ + H^^2^^) / H^^^\ is found to be 1.15, which also matches well with its value reported 
by Metlushko et al. [2^ . 

In Figs. 5 (a) and 5 (b), the field region is shown to be subdivided into three parts, (i) H < Hon 
, (ii) Hon < H < Hp and H > Hp, and designated as well ordered Bragg glass (BG) state, 
polycrystalline vortex glass state and amorphous state, on the basis of notions articulated in 
the literature [21 HI [30]. The anisotropy in the phase diagram in Fig. 5 is striking, and may 
appear perplexing at the first glance. We rationalize a subtle difference between those two 
phase diagrams by recalling that the spatial order prior to Hon for H || [110] is superior to that 
for H [| [100]. Thus, on progressively tilting the field from [110] towards [001], the field region 
of BG state gradually shrinks and reaches a minimum value at 6 = 60° (cf. Fig. 5(a)). On 
the other hand, starting from [100], the minimum is reached earlier at 6 = 30° (cf. Fig. 5(b)) 
and, thereafter, the field interval of BG enhances as Hon values climb up towards its plateau 
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value for H || [001]. Between Hon and Hp, the vortex state is expected to comprise coexistence 
|32j of ordered and disordered phases and above Hp, the disordered phase is the underlying 
equilibrium state, as had been shown for experiments in V^Si |26j and 2H-A^6S'e2 [32j. In the 
field ramping mode of experimentation, it has been shown [331 [34] that disordered bundles of 
vortices injecting through edges and corners continuously attempt to anneal to the underlying 
equilibrium vortex state, thereby, implying that residual disorder most likely coexists with the 
ordered state even prior to Hon- In platelet shaped sample, like the LNBC crystal, the extent of 
disorder injecting into the sample through edges and surfaces could depend on the orientation 
of the crystal, as has been witnessed while studying SMP phenomenon in crystals of high Tc 
Bismuth cuprate, which had been cut in the form of a platelet and a prism [35]. While such a 
possibility can be a source of anisotropy in disorder /pinning in our case as well, we are more 
inclined to believe in the correlation between the anisotropy in disorder and the anisotropy in 
vortex lattice symmetry transition. The anisotropy in pinning observed between [110] and [100] 
in LNBC is analogous to that observed in a spherical crystal of (cubic) V^Si [26]. The shape 
dependent argument would not suffice to rationalize observations in spherical sample of V^Si. 



Conclusion 

The observation of large angular variation in the vortex phase diagram of a superconducting 
system, which prime facie imbibes marginal anisotropy in its intrinsic superconducting param- 
eters, is a curious finding. Usually vortices are expected to assume triangular symmetry (apex 
angle /? = 60°) in the mixed state of a type-II superconductor, however, in quaternary boro- 
carbide superconductors, YNi2B2C and LuNi2B2C, the local symmetry of vortex lattice has 
been known to change from rhombohedral (/? < 60°) towards square {P = 90°) at low fields. 
We note that for vortex arrays generated by applying field in different crystallographic direc- 
tions of tetragonal structured (Y/L)NBC, the field interval over which good spatial order in 
the vortex lattice prevails correlates with the details of the underlying VL symmetry transition. 
For field oriented along crystal directions of higher symmetry, like, the c-axis or the ab-basal 
plane, where VL symmetry becomes perfect square or quasi-square, the well ordered Bragg 
glass state transforms to the disordered phase of vortex matter at high fields. Also, the c-axis, 
along which the VL assumes square symmetry, the correlation volume for the spatially ordered 
state is the largest and the shrinkage in the correlation volume across order-disorder transition 
a la peak effect is maximum. When the applied field is oriented away from c-axis towards the 
ab-basal plane, at intermediate angles {6 ~ 30° to 60°), the field interval of Bragg glass state 
shrinks presumably due to enhancement in effective pinning. The fact that application of a 
large driving force can overcome the enhanced effective pinning and restore the field interval 
of Bragg glass phase supports the above assertion. Detailed field-angle dependent studies in 
suitably shaped (L/Y)NBC crystals are desired to fortify the correlation between anisotropy in 
VL symmetry transition and the angular variation in their vortex phase diagram. Experimental 
results presented will hopefully spur theoretical studies to quantitatively account for them. 
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